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Abstract 

Operational performance, algorithms, stability and physics results of the Outer em calorime- 
ter of FOCUS are overviewed. 
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1 Introduction and historical overview 

FOCUS is a heavy-flavor photoproduction experiment located at the Wide Band Photon 
Area of Fermilab. An upgraded version of its predecessor E687[ p, FOCUS is composed 
of (Figjj]) a strip silicon detector, a large acceptance magnetic spectrometer with MWPC, 
Cherenkov differential counters, muon detectors, hadron calorimetry, and both forward 
(Inner em), and large angle (Outer em) electromagnetic calorimetry. The topics in charm 
physics that are being investigated include lifetimes, semileptonic decays, charm baryons, 
charm spectroscopy, searches for D° — D° mixing and for rare and forbidden decays, 
charmonium production and radiative decays, and charm meson and baryon decays with 
neutrals. 

The experiment was designed starting in 1981; in 1985 there was the first test beam 
of the completed spectrometer; in 1987 the first E687 data-taking period, interrupted by 
a fire. The spectrometer was seriously damaged, but quickly repaired. In 1990-1991 the 
second E687 data taking took place. In 1995 the E687 spectrometer became FOCUS and 
underwent upgrades, with the Outer em calorimeter being equipped with a new plane of 
square scintillator tiles. The 1996-1997 FOCUS data-taking period met the goal of col- 
lecting ten times the E687 statistics of reconstructed charm decays, by fully reconstructing 
more than ten million charm particle decays. 

2 Physics requirements, geometry acceptances 

The Outer em calorimeter (Fig. 0a) is located 900 cm from the target. Its external di- 
mensions are (255 x 205) cm 2 , with an internal rectangular aperture (51 x 88) cm 2 . This 
corresponds to an angular acceptance for photons (28 < \9 X \ < 142) mrad, (49 < \9 y \ < 
114) mrad. A vertical gap, set at 9 cm for the 1996-97 run, avoids showers from the most 
abundant background process, i.e., Bethe-Heitler e + e pair production. The Outer em 
is required to reconstruct 7-initiated showers from charm radiative and tt° decays in the 
energy range (0.5 < E 1 < 15) GeV and to perform e/n identification for charm semi- 
electronic decays in the momentum range (2.5 < P < 20) GeV/c, thus extending the 
Cherenkov counter identification, which is limited to P < 6 GeV/c. Some fi/n identi- 
fication power is expected to help identify low-momentum muons in charm semimuonic 
decays. Typical geometrical acceptances range from 30% for electrons and muons in 
charm meson semileptonic decays, and 40-50% for decays with one or more 7r° in the 
final state, including the case of shared n° (one 7 in the Inner and one in the Outer em 
calorimeter). 
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Figure 1 : The FOCUS (E83 1) spectrometer at the Fermilab Wide Band Photon Beam. 
3 Mechanical structure 

The Outer em calorimeter is made of Pb plates (stiffened with a 6% Sb by weight) and 
scintillator layers (POPOP, C24H42N2O2 doped with 8% naphtalene, and NE-102 were 
used), for a total of 19 X and 1.4 Aj (Tab.[TJ). Scintillator layers are made of strips, whose 
light readout is either individual (OEO, OE9 segments) or five-fold integrated by a light 
guide to a single PM (all other segments), for a total of 1136 readout channels. Each 
counter in the calorimeter is individually wrapped in 0.1 mm Al foils and black plastic. 
Horizontal and vertical five-fold counters are interlaced as shown in FigfJ b). The coun- 
ters are arranged in nine independent views along Z (Z is the beam direction), and four 
independent quadrants in the (X,Y) plane. A module of S-Z strips (45° — 135°) performs 
horizontal-vertical matching of clusters. The counters are equipped with ten-stage, EMI- 
9902KB photomultiplier tubes (PMT) operating at a typical gain of 10 6 at 1000 V, with 
a quantum efficiency of 20% at 440 nm, which were individually tested in order to select 
only those with good linearity and small sensitivity to rate effect[ Q|. The PMT's are pow- 
ered by LeCroy 1440 and custom-made FRAMM[ |3|] HV systems, via a high-linearity, 
anode grounded voltage divider supplying 1.5 mA at 1500 V. The PMT signals reach 
the counting room via 60 m long coaxial cables, where they are converted by a Lecroy 
1881M Fastbus ADC, with a 0.050 pC/count conversion. The Outer em calorimeter can 
be displaced both horizontally and vertically for calibration and access purposes. 

A scintillator tile array module recovers showers in the small-angle, high-occupancy 
region, improves horizontal-vertical matching, and cleans the n° peak by rejecting fake 
matches. The module is located at shower max, i.e., between the OE1 and OE2 modules, 
and is composed of an array of 100 supertiles, and edge counters to flag laterally noncon- 
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Figure 2: a) The Outer em calorimeter longitudinal structure showing the sandwich of 
absorber and scintillator planes, b) X and Y strip planes interlaced and five-fold integrated 
to a readout PMT (dimensions are millimeters). 



4 



Table 1 : Longitudinal segmentation and counter geometry of the Outer em calorimeter. 
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Figure 4: left) A complete tile counter before wrapping; right) layout of tile counter 
array, and edge counters. 



tained showers. Each supertile is made of two (10x10x0.5) cm 3 tiles (BICRON BC404- 
B scintillator), each is equipped with two 20 cm-long, 1-mm-diameter wavelength shifting 
(WLS) optical fibers (Kuraray Yl 1 multiclad S-type) 20 cm-long, a-cut, with all four ends 
thermally spliced to 2 m-long clear fibers following the CDF endplug splicing technique 
[ Q], with a heat- shrinking tube to protect the splice. The eight ends of the clear fibers of 
each supertile are optically coupled by means of optical grease to a EMI-9902KB PMT. 
Each tile is wrapped in white Teflon tape, and the side is 
painted with white reflective paint by BICRON. The tile 
array is enclosed in a light-tight Al case. The light trans- 
mission efficiency of the thermal splice was measured on 
relevant samples during the splicing process and shown 
to be typically 94%. 
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4 Equalization with mips 

Beam halo muons were used to determine the counter ge- 
ometry, for an approximate PMT gain balancing (±10%), 
to determine the light attenuation curves inside strips 
(Fig.|]), and the equalization constants. The light output 
was measured in the laboratory with cosmic rays to be 
30 photoelectrons/mip for a supertile, and 100 photoelec- 
trons/mip for a five-fold counter at 20 cm from PMT. The 
relative width fwhm of the energy distribution for a mip 




Figure 3: Light attenua- 
tion curves for three different 
sized strip counters. 
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Figure 5: (a) Loss of reconstruction efficiency in the small-angle region of the Outer em 
calorimeter; (b) 77 invariant mass peak with and without tile array hit confirmation. 

is 40% for a supertile, and 30% for a five-fold counter. 

5 Shower reconstruction strategy 

The reconstruction algorithm begins with the identification of clustered energy deposits. 
Clusters of energy deposits associated with the projection of charged tracks reconstructed 
in the magnetic spectrometer are tagged. The reconstruction of neutral showers uses the 
remaining clusters. Energy deposits (in ADC counts units) for each counter are multiplied 
by the mip equalization constants, and then summed up to determine the detected energy 
associated to each cluster. 

Pairs of clusters in the two orthogonal X-Y views of each OE segment are formed 
using the energy balance as a criterium, after proper weighing for light attenuation inside 
the strip counters. Neutral showers are formed by aligned X-Y pairs in the different OE 
segments. The diagonal counters and the tiles are used to resolve ambiguities. The use of 
tile array information improves the efficiency of X-Y matching in the small angle region 
(Fig.|5]a) by reducing the number of fake matches (FigJJb). 

The coordinates of the shower centroids are determined considering the energy de- 
posited in each counter of the X-Y pairs. Once corrected for systematic effects (§ 6), 
the shower centroids determine the photon incidence point. The n° invariant mass is 
computed by the measurement of relative angle and energy of the decay photons. When 
computing the invariant mass of higher states with one or more n° in the final state, the n° 
invariant mass is fixed at its nominal rest value, and the photon momenta are rearranged 
by means of a 1-C fit. The classical algorithm in Ref.[ [7|] was modified in order to take 
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into account space resolution[|8p. 

Neutral showers are identified as em or hadronic by means of the Discriminant 
Analysis algorithm (§ 7). Finally, the sum of energy clusters longitudinally forming a 
reconstructed photon track gives the detected energy E detected . The reconstruction effi- 
ciency for single isolated showers was measured using primary Bethe-Heitler e + e~ pairs. 
The e + e tracks found in the proportional chambers were projected onto the OE front. 
Shower reconstruction was then performed using all the available clusters. The efficiency 
for reconstructing the shower associated with the electron or positron track was greater 
than 95% over the range 2-20 GeV, and better than 90% over the range 0.5-2 GeV. 

6 Energy calibration, linearity, energy and space resolution 

The response of the calorimeter to photon- and electron-initiated e.m. showers, the scale 
factor a between detected energy and incident energy, and the energy resolution have 
been studied using a GEANT simulation. Simulation predictions have been verified with 
e + e~ pairs and electron beam in calibration runs, e + e~ pairs from the process n° — > 
77, jj\f — > e + e~ of photon conversion in physics events for electron-initiated showers 
constants , and n° peaks for photon-initiated showers constants and absolute calibration 
of the energy scale. 

The detected energy is parametrized as Edetected = E/a where E is the particle 
incident energy and E det ected is the particle energy deposited in the calorimeter active lay- 
ers. Energy linearity and resolution are shown in Figj£], in agreement with the simulation 
predictions. 

The photon impact point X true is determined from the em shower center-of-gravity 
X cog = 2AJ2 iXiAi/ iAi, where 2A = 3.3 cm is the counter width, after applying the 
standard correction[ ||] X true = frarcsinh (^psinhy) After correction, we determine 
the space resolution as from e + e~ calibration events (Fig.0) to be a(X true ) = ±0.3 cm in 
the energy range (1 < E 1 < 20) GeV. 

7 Particle ID 

The Outer em calorimeter extends the e/n rejection beyond the Cherenkov momentum 
range, i.e., from 6 to 20 GeV/c. The identification algorithms have been developed in the 
framework of the Discriminant Analysis [ §], which allows one to distinguish between two 
or more groups of events. As first step, we determine a set of N variables V\ : n (Discrimi- 
nant Variables), significantly different among the M groups of events ({A,}, j = 1, M) 
to be distinguished (Fig||). A typical set of variables is composed of the ratio E/P be- 
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Figure 6: Energy resolution (with the constant term added in quadrature) and linearity as 
measured by calibration e + e~ pairs. 

tween the energy measured by the calorimeter and the track momentum by the track- 
ing system, the lateral and longitudinal shower development pattern, the cluster centroid 
residuals and widths. Next, the score function Sa = Z)i=i n ls built, and we find coef- 
ficients q which maximize separation among the Sa, thus applying one cut on S (Fig.|9|a). 
As training samples of known membership we used e + e~ Bethe-Heitler pairs embedded 
in hadronic events, pions from K s — > n + n~ decays, and muons from dedicated runs with 
beamdump. The overall pion residual contamination obtained for an 85% electron effi- 
ciency is 10~ 2 , while the pion residual contamination is 10 _1 for 85% muon efficiency 
(Fig.Pb). Efficiency for muons from J/ip decay, and rejection of pions compared with the 



Outer muon detector performances are shown in Fig.l 1 a,b). 
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Figure 7: a,b) Shower center-of-gravity X cog as measured for a beam of calibration 
electrons sweeping the Outer em at increasing angles. The calorimeter response degrades 
when the em shower reaches the lateral boundary of the detector acceptance, thus losing 
full lateral containment, c) Distribution of residuals after correction. 

8 Long-term stability 

The calorimeter stability is controlled by monitoring the ADC pedestals between spills, 
the PMT HV supplies, and the single-channel rough response by a N 2 laser light source[ 
|9j]. The availability of a muon beam halo over the entire area of the Outer em allows 
muon calibrations to be performed regularly. The long-term, fine-grained run-dependent 
stability is given by exploiting physical signals in events, namely 

1. 7T° -> 77 (Fig.[lO| a) 

2. E/P for electrons (e + e~ Bethe-Heitler pairs embedded in hadronic events, Fig.[T0|b) 

3. 7T° — > 77, with one 7 conversion 7JV — > e + e~ (Fig.[To|c) 



Results are summarized in Fig.[lO|d,e,f. Electron and n° signals in hadronic events 
(Fig.[lO|b,a) can track the shifts in detector response up to a stability of ±1% over the entire 
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18-month data-taking period. Work is in progress to study process (3) (Fig.|TTS| c), which 
will provide on an event-by-event basis the photon energy calibration constant, electron 
and positron E/P, and a check of the stability of the measurement of P, as performed by 
the magnetic spectrometer. 

9 Physics results 

Clean signals of the decays of charm mesons and baryons in many channels have been 
obtained by using the Outer em electron and muon identification and the n° reconstruction 
(Fig®. 

The geometrical acceptance of the Outer em for electrons and muons from semilep- 
tonic charm decays is about 30%, extending especially in the low-momentum region 
P ~ 4 GeV/c. The detection of charm semileptonic decays with the electron momen- 
tum in the acceptance of the Outer em is of particular interest in view of the possibility 
of determining the q 2 dependence for the formfactors (see, eg, Ref.[ [RJ for an updated 
review). Studies performed in E687 include the Cabibbo-suppressed semileptonic decay 
D° — > 7r - e + z/ e (Ref.[ [TT]]). The efficiency of fi/n rejection was measured on J ftp de- 
cays where the muons were identified by the Outer muon counters, and found to be better 
than 80% from 5 to 50 GeV/c (Fig.|TT|a). For the fi/n rejection, similarly to the case of 
e/ir, the contribution of the Outer em is effective especially in the low momentum region 
(Fig{TT]b), as measured on pions from K s decays. 

Numerous charm decays have been found with 7r°'s reconstructed by the Outer em 
(Fig.pT]c-m). Thanks to the large statistics and low level of background, precision mea- 
surements such as the isospin mass splittings m(D* + ) — m(D*°) of excited charm meson 
states have been initiated. Preliminary results[ [T3p (Fig.|lT)s) show how well the attain- 
able precision compares with the best results obtained by e + e~ experiments using crystal 
calorimetry[|l4]]. 

In conclusion, the ten-year operational experience of the Outer em calorimeter of 
FOCUS shows how a detector based on conventional techniques is able to perform con- 
sistently and provide competitive physics results. The implementation of a scintillator 
tile tiebreaker has increased shower reconstruction efficiency in the small-angle region, 
and has considerably cleaned the ir° peak of spurious combinations. Long-term response 
stability of 1% is attained by cross-calibration between E/P for electrons and n° peaks in 
hadronic events, while Discriminant Analysis is used to provide e/n and fi/ir identifica- 
tion. 
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Figure 8: Some Discriminant Variables for e, n and ji. 
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Figure 9: a) Distribution of the score variable S for electron, muon, and pion train- 
ing samples; b) efficiency and residual contamination for electron/pion and muon/pion 
identification. 
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Figure 10: Long-term detector stability over 18 months, (top, left to right) 77 invariant 
mass distribution with the n° peak, E/P distribution, and 7e + e~ invariant mass distribution 
with the 7T° peak, (bottom) Timelines, as a function of run number (typically one run 
per hour), of the E/P and tt° peak, and the E/P peak corrected by the it peak. After 
correction, the detector response is stable within ±1%. The residual shift is compatible 
with the stability of the charged track spectrometer. 
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Figure 11: FOCUS preliminary results on a fraction of the 1996-97 dataset. From 
top to bottom: efficiency on muons from J/ip — > ji + n~ and rejection of pions from 
K s — > 7r + 7r~; D*° and D* + decays with 7r° in the final state, and comparison of the pre- 
cision attainable on the isospin mass splittings compared to world average, as functions 
of detached vertex significance I /a; several decays of D + and D° mesons with ir° in the 
final state, with a selection of preliminary background-reducing cuts. 
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